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This paper reports on the most comprehensive data set obtained on differential and fully integrated 
cross sections for the process ep — » e'pir + ir~ . The data were collected with the CLAS detector at 
Jefferson Laboratory. Measurements were carried out in the so-far unexplored kinematic region of 
photon virtuality 0.2 < Q 2 < 0.6 GeV 2 and invariant mass of the final hadron system W from 1.3 
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states in kinematical dependencies of the 1-fold differential cross-sections. 
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I. INTRODUCTION 

An extensive research program is currently underway 
in Hall B at Jefferson Lab with the CLAS detector, fo- 
cused on studies of nucleon resonances (N*) in various 
exclusive channels of meson electroproduction off pro- 
tons P, H, 0]. Part of this effort is aimed at the de- 
termination of electrocouplings for an entire spectrum 
of excited nucleon N* and A* states versus the photon 
virtuality Q 2 = — (e — e') 2 , where e and e' are the 4- 
momcntum vectors of the incoming and scattered elec- 
tron, respectively Comprehensive information on res- 
onance electrocouplings and their evolution with Q 2 is 
needed to probe the spatial and spin structure of the 
resonance transitions. This information is needed to en- 
hance our understanding of the effective strong interac- 
tion that is at the core of internal baryon structure and 
the decay of baryons. It is also needed to firmly estab- 
lish the connection of effective degrees of freedom such 
as a) 'constituent' quarks in the binding potential or b) 
'constituent' quark scattering amplitudes to the elemen- 
tary quarks and gauge gluons of QCD, the theory of the 
strong interaction. 

In the past five years, single pseudoscalar meson elec- 
troproduction has been studied in several exclusive pro- 
cesses, e.g. pn°, nir + , prj, KA, and KT,, 0, H, @, 0, 
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1, 1, M EL El M M HE HEEim, which included 
differential cross sections with complete polar angle and 
azimuthal angle distributions, as well as several polariza- 
tion observables. From these data sets, resonance transi- 
tion electromagnetic form factors have been determined 
for the A(1232)P 33 , AT(1535)5ii, and iV(1440)Pii, cov- 
ering a wide range in Q 2 . 

Studies of charged double pion electroproduction off 
protons represent an important part at this effort. Sin- 
gle and double pion production are the two largest con- 
tributors to the total photo- and electroproduction cross 
sections off protons in the resonance region. The final 
states produced in these two exclusive channels have con- 
siderable hadronic interactions, or so-called final state in- 
teractions (FSI). FSI may be determined using the data 
of experiments with hadronic probes [l9l |. According to 
these data, the cross section for the irN — ► 7r7rA^ reaction 
has the second largest strength of all of the exclusive 
channels in the ttN interaction. Considerable FSI be- 
tween the irN and irirN final states result in substantial 
contributions to the amplitudes of both single and double 
pion electroproduction from the electroproduction ampli- 
tudes of the other channel. Accounting for these coupled- 
channel effects is essential in order to get the amplitude 
description compatible with the constraints imposed by 
unitarity. Therefore, for N* studies both in single and 
double pion electroproduction, information is needed on 
the mechanisms contributing to each of these channels 
in order to take properly into account the impact from 
coupled-channel effects on the exclusive channel cross sec- 
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tions. The knowledge of single and double pion electro- 
production mechanisms becomes even more important 
for N* studies in channels with smaller cross sections 
such as prj or KA and KY, production, as they could be 
significantly affected in leading order by coupled-channel 
effects produced by their hadronic interactions with the 
dominant single and double pion elcctroproduction chan- 
nels. Therefore, comprehensive studies of single and dou- 
ble pion electroproduction are of key importance for the 
entire N* research program. 

The world data on double pion electroproduction in the 
nucleon resonance excitation region were rather scarce 
before the data from CLAS became available. Fully in- 
tegrated cross sections for it A isobar channels as a func- 
tion of the invariant mass of the hadronic system W 
and photon virtualities Q 2 were available from DESY 
[20| . However, the data are presented in very large kine- 
matic bins, AW = 200 - 300 MeV and AQ 2 = 0.2 - 0.6 
GeV 2 . Center-of-mass angular distributions for ir~ were 
also measured, but were averaged over a very large inter- 
vals in Q 2 from 0.3 to 1.4 GeV 2 . This makes it virtually 
impossible to determine the nucleon resonance parame- 
ters from such measurements. 

The first detailed data on charged double pion elec- 
troproduction cross sections in the resonance region were 
obtained with CLAS [U The data were collected 
for W=1A - 2.1 GeV in 25 MeV bins and for Q 2 = 
0.5 — 1.5 GeV 2 in 0.3 GeV 2 wide bins. The current exper- 
iment covers invariant masses of tt~tt + , ir + p and ir~p in 
each (W, Q 2 ) bin for Q 2 = 0.2 - 0.6 GeV 2 in 0.05 GeV 2 
wide bins and for W =1.30 to 1.57 GeV with 25 MeV 
bins. In addition, angular distributions for 7r~, 7r + , and 
proton, as well as angular distributions in q, (i=l,2,3) 
angles (see Sect. |V]for on definitions), were measured. 
These very detailed measurements are crucial to deter- 
mine the most significant production mechanisms for this 
process. 



II. ANALYSIS TOOLS 

The presence of three hadrons in the final state 
presents considerable complications in the phenomeno- 
logical analysis. Efforts to apply partial wave analysis 
(PWA) techniques to double pion production by electro- 
magnetic probes are limited to photoproduction, where 
very high statistics data are available [25|, HH, [27[. A 
strong reduction in statistics for individual bins in Q 2 
makes application of PWA methods in double pion elec- 
troproduction data much more difficult. Moreover, there 
is no model-independent way to disentangle resonant 
and non-resonant mechanisms in any given partial wave. 
Therefore, reaction models are needed to isolate the reso- 
nant parts in the double pion production amplitudes and 
evaluate the N* electromagnetic transition form factors. 

Following the pioneering effort of Ref . [28| , several ap- 
proaches have been developed more recently for the de- 
scription of double pion photo- and electroproduction in 



the resonance region [H, [H [M H H [H, [H[ . These 
efforts were based on a very limited amount of experi- 
mental data: mostly on W and Q 2 dependencies for fully 
integrated cross sections and on invariant mass distribu- 
tions for various pairs of the final state hadrons. The re- 
action models used mcson-baryon degrees of freedom: N, 
A, 7r, a and p. Effective meson-baryon Lagrangian opera- 
tors were constructed based on Lorentz invariance, gauge 
invariance and crossing symmetry. For the description of 
experimental data, a limited set of non-resonant meson- 
baryon diagrams was used with amplitudes calculated 
from effective Lagrangians together with contributions 
from several, mostly low-lying nucleon resonances (M < 
1.6 GeV). A general framework for the implementation 
of other meson-baryon degrees of freedom was proposed 
in Ref. [HJ , but so- far has not been fully realized. 

The meson-baryon diagrams in reaction models may 
account for many partial waves. However, in any reac- 
tion model we need to truncate the infinite set of meson- 
baryon diagrams, keeping just the relevant mechanisms. 
Moreover, the choice of a particular effective Lagrangian, 
describing meson-baryon interactions, may be done only 
at a phenomenological level. 

At the distance scale appropriate for the size of 
hadrons, the amplitudes of effective meson-baryon inter- 
actions contributing to the reaction cannot be expanded 
over a small parameter except for a small kinematic re- 
gion near threshold W and Q 2 < 0.2 GeV 2 accessible 
for chiral perturbation theory. This feature makes it im- 
possible to select contributing diagrams based on a per- 
turbative expansion for the entire double pion reaction 
phase space covered by the CLAS measurements. So-far, 
no approach has been developed based on a fundamental 
theory that would allow either a description of an effec- 
tive Lagrangian or a selection of the contributing meson- 
baryon mechanisms from basic principles. We therefore 
have to rely on fits to the now available detailed experi- 
mental data sets to develop reaction models that contain 
the relevant mechanisms. 

The large reaction phase space coverage of CLAS data 
opens up qualitatively new opportunities for the analy- 
sis of charged double pion electroproduction. The ex- 
clusive channel ep — * e'pTt + -K~ offers many observables 
for the analysis. The hadronic final state can be pro- 
jected on nine independent 1-fold differential cross sec- 
tions in each W and Q 2 bin. For the first time, all of 
these observables have become experimentally accessi- 
ble [U, ■ By studying the kinematical dependencies 
of the differential cross section and their correlations we 
are able to establish the presence and strength of the 
relevant reaction mechanisms. A phenomenological ap- 
proach 0, m m, m m no, m s, ei w as developed 

in collaboration between Jefferson Lab and Moscow State 
University, referred to herein as " JM" . This approach is 
intended to establish all significant mechanisms seen in 
the observables of charged double pion electroproduction, 
to isolate the resonant parts of the amplitudes, and to de- 
rive the electrocouplings of nucleon resonance transitions 
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from fits of all measured observables combined. 

The CLAS data presented in this paper were collected 
with binning resolutions over W and Q 2 surpassing by 
almost an order of magnitude what was achieved in pre- 
vious measurements before the experiments with CLAS. 
These data will allow us to extract the electrocoupling 
amplitudes of the A(1440)Pn and A(1520)L>i 3 states. 
A certain advantage of the double pion channel is that 
the amplitudes of the iV(1440)Pu resonance do not in- 
terfere with the high-mass tail of the A(1232)p33 state, 
as is the case for the amplitudes in single tt production. 

The behavior of N* clectrocouplings at small photon 
virtualities is of particular interest. Studies of the mag- 
netic transition form factor for the A(1232)p33 in Ref.|45j 
revealed considerable meson-baryon dressing effects in 
addition to the 3-quark core contributions. The dress- 
ing is expected to decrease with increasing Q 2 . It is 
most pronounced at Q 2 < 1.0 GeV 2 Hl|. While the 
role of meson-baryon dressing effects has been studied 
based on the data on electrocouplings of the A(1232)p53 
state, for resonances heavier than the A(1232)P 3 3, this 
remains an open question, and is currently being ad- 
dressed at Jefferson Lab through extensive theoretical ef- 
forts [46j, |47|, |48j, |49| . Accounting for these effects is a nec- 
essary step to probe quark and possibly gluonic degrees of 
freedom in baryons. The comparison of constituent quark 
model predictions [EH, [E3 with the measured N* elec- 
trocouplings [1, Q , as well as the coupled-channels anal- 
ysis of the data on reactions with hadronic probes [49| . 
suggest considerable meson-baryon dressing effects at Q 2 
< 0.5 GeV 2 for the A(1440)Pn and A(1520)P>i 3 elec- 
trocouplings. Therefore, the information on these states 
at low photon virtualities from charged double pion elec- 
troproduction may further elucidate the relevant degrees 
of freedom in resonance excitation at hadronic distance 
scales. 

The data presented in this paper can also provide in- 
formation on the p — > A axial transition form factor. 
Current algebra [53| relates the contact term in the set 
of non-resonant Born terms for tt A isobar channels to the 
axial transition form factors. These contact terms could 
be fit to the data within the framework of the JM ap- 
proach [4l|, El| . So-far, the p — » A axial transition form 
factor has been determined mostly from neutrino-induced 
reactions f54l. l55l. [56j . Axial transition form factors offer a 
complementary view of baryon structure, seen in the ax- 
ial vector currents, while electroproduction experiments 
usually access baryon structure through vector currents. 
Recently, lattice QCD results have become available [57| , 
which make the experimental study of the nucleon axial 
structure an important topic of hadronic physics. 



III. EXPERIMENT 

The measurement was carried out using the CE- 
BAF Large Acceptance Spectrometer (CLAS) [Hj] at the 
Thomas Jefferson National Accelerator Facility (Jeffer- 



son Lab). CLAS provides almost complete angular cov- 
erage in the center-of-mass frame. It is well suited for 
conducting experiments that require detection of two or 
more particles in the final state. Such a detector and the 
continuous beam produced by CEBAF provide excellent 
conditions for measuring the ep — > e'p'ir + ir~ cross sec- 
tion by detecting the outgoing electron, proton and at 
least one pion in coincidence. 



A. Apparatus 

The main magnetic field of CLAS is provided by six su- 
perconducting coils, symmetrically arranged around the 
beam line, which generate an approximately toroidal field 
in the azimuthal direction around the beam axis. The 
gaps between the coil cryostat are instrumented with six 
identical detector packages, referred to here as "sectors", 
as shown in Fig. [TJ Each sector consists of three regions 
(Region 1, Region 2, Region 3) of Drift Chambers (DC) 
[59j to determine the trajectories of the charged parti- 
cles as they travel from the target outward in magnetic 
field, a Chcrenkov Counter 



6C | for electron identification, 
Scintillator Counters (SC) [6l|] for charged particle iden- 
tification using the Time-of-Flight (TOF) method, and 
an Electromagnetic Calorimeter [62] for electron identi- 
fication. The liquid-hydrogen target was located in the 
center of the detector. To reduce the electromagnetic 
background resulting from M0ller scattering off atomic 
electrons, a second smaller normal-conducting toroidal 
magnet (mini-torus) was placed symmetrically around 
the target. This additional magnetic field prevented 
M0llcr electrons from reaching the sensitive detector vol- 
ume. A totally absorbing Faraday cup, located at the 
very end of beam line, was used to determine the in- 
tegrated beam charge passing through the target. The 
CLAS detector provides k 80 % of 4-7T solid-angle cov- 
erage. The efficiency of detection and reconstruction for 
single stable charged particles in the fiducial regions of 
CLAS is greater than 95 %. The combined information 
from tracking in the DC and the SC systems allows us to 
reliably separate protons from positive pions. Additional 
constraints for event selection come from the overdeter- 
mined kinematics, which allows use of the missing mass 
technique. 

Due to possible slight misalignments in the DC po- 
sitions and small inaccuracies in the description of the 
torus magnetic field, the reconstructed momentum and 
angle of particles may have small systematic deviations 
from the physical value. To correct these deviations, elas- 
tic electron-proton scattering was checked and the elec- 
tron 3-momenta were corrected to ensure the proper mass 
peak position for the recoil proton [63| . The proton en- 
ergy losses in CLAS were estimated from a simulation 
of proton propagation through the detector materials in 
kinematics corresponding to charged double pion electro- 
production. 
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FIG. 1: (color online) Cross-sectional views of the CLAS de- 
tector. The top panel shows a cut along the beam line and 
through the mid-plane of two opposite sectors. The bottom 
panel shows a cut perpendicular to the beam line and through 
the nominal target center. Descriptions of the detector ele- 
ments are given in the text of Sect. IIII Al 



B. Data taking and data reduction 

This analysis is based on data taken during the 1999 
elc run period. The 1.515 GeV electron beam at a cur- 
rent of 3 nA was incident on a 5-cm-long liquid-hydrogen 
target corresponding to an instantaneous luminosity of 
~ 4 x 10 33 cm~ 2 s~ 1 . The size of the beam spot at 
the target was ~ 0.2 mm, with position fluctuations < 
±0.04 mm. The main torus current was set at 1500 A, 
which created a magnetic field of about 0.8 T at polar 
angles of 20° that decreased with increasing polar an- 



gle. The CLAS event readout was triggered by a co- 
incidence of signals from an electromagnetic calorimeter 
module and a threshold gas Cherenkov counter in one of 
the six sectors, generating a total event rate of ~ 2 kHz. 
The number of accumulated triggers at these detector 
settings was about 4.2 x 10 s . These data were further 
analyzed to extract the differential cross sections for the 
ep — > e'p'Tt + TT~ reaction. 



IV. EVENT SELECTION 

The ep — > e'p'TT + Tr~ reaction is selected by measuring 
the scattered electron, as well as the proton and ir + in the 
hadronic final state. In the magnetic field configuration 
used in this measurement, the negatively charged pions 
have a smaller probability for detection than the posi- 
tively charged particles. However, the process is kine- 
matically overconstrained and the detection of all parti- 
cles in the final state is not required for an unambiguous 
identification of the exclusive reaction. To retain max- 
imum acceptance, we chose to not require detection of 
the n~ , but infer the presence and kinematics of the un- 
detected 7T~ by computing the mass of the undetected 
particle from 4-momentum conservation and its charge 
from charge conservation. The two other topologies, with 
either undetected proton or undetected ir + , have consid- 
erably lower rates and were used for systematics studies 
and for cross checks, but arc not included in the deter- 
mination of the final cross sections. 
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FIG. 2: (color online) Distribution of photoelectrons in the 
Cherenkov counter for 0.3 < Q 2 < 0.4 GeV 2 in two W bins 
centered at the values shown in the plot. The curves represent 
a Poisson fit. 



For each event we identify as the electron candidate the 
first coming negatively charged particle detected in the 
electromagnetic calorimeter and the Cherenkov counter. 
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To select true electrons, we apply a cut in the number of 
photoelectrons (N p h e > 2.5) produced by the Cherenkov 
light signal in the photomultiplicrs. This cut also elimi- 
nates a small fraction of electrons (< 6 %), as shown in 
Fig. [2 The shadowed areas correspond to the cut-out 
electrons. A special procedure was developed to account 
for these electrons in the evaluation of the reconstruction 
efficiency, based on the extrapolation of the photoclcc- 
tron spectra into the cut-out areas using a fit based on a 
Poisson distribution. The quality of electron identifica- 
tion may be seen in Fig.[3j where we display the energies 
deposited in the outer part of the calorimeter versus the 
energies deposited in the inner part of the EC, normal- 
ized to the momenta of the incoming particles. A spot 
from minimum- ionizing pions, clearly seen in inclusive 
electron events (top part of Fig. [3]) disappears after ap- 
plying the described cuts for electron selection (bottom 
part of Fig. El). 



for positively charged particles. 
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FIG. 3: (color online) Distributions for the energies deposited 
in the outer part (Eout) versus the energies deposited in the 
inner part (Ei n ) of the EC normalized to the momenta of the 
outgoing particles (P e ')- The distribution for the accumulated 
triggers is shown in the top plot, while the distribution for the 
events selected, applying the photoelectron cut, is shown in 
the bottom plot. 



Using information from the time-of-flight scintillators 
and the path length determined by tracking, the parti- 
cle's velocity (/3) was determined. The information from 
the drift chambers, combined with the known magnetic 
field, provides a measurement of the particle momentum 
(p). Relativistic relations between particle mass, momen- 
tum and velocity were used in order to determine the 
particle's mass, allowing us to separate pions, kaons, and 
protons in the kinematic range covered by the measure- 
ment. Fig. 2] shows the velocity versus the momentum 
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FIG. 4: (color online) j3 versus momentum for positively 
charged hadrons; solid curves show the cuts used to identify 
pions and protons; dashed curves show the explicit j3 vs p re- 
lationships, using the pion and proton mass, respectively. The 
events near (3=1 and p=0.1 GeV are positrons, e.g. from ty° 
Dalitz decay. The data were collected in a single scintillator 
bar. 



After identification of the three particles, events were 
selected with one electron, one proton and one n + . The 
squared missing mass distribution is shown in Fig. 
which clearly shows the pion mass peak and also indi- 
cates that multi-pion background (> 2tt) contributes less 
than 1 % to the total number of charged double pion 
events selected by the cuts. The small background is re- 
lated to the kinematic coverage of our experiment W < 
1.57 GeV, where the exclusive channels with more than 
two pions in the final state are suppressed due to their 
thresholds. The almost negligible contribution from the 
multi-pion background and the rather small radiative ef- 
fects (see Sect. IV F| allowed us to apply a wide exclusivity 
cut over the squared missing mass distribution, shown in 
Fig. [5] by the two arrows, in order to collect the majority 
of the 2tt events. 

The CLAS detector has an active detection solid angle 
smaller than 4n due to the space filled with the torus 
magnet coils. The angles covered by the torus magnet 
coils are not equipped with any detection system and 
therefore give rise to inactive areas. The boundaries of 
the active areas are not well defined and do not provide 
regions for particle reconstruction with full reconstruc- 
tion efficiency. Therefore, for the analysis we accept only 
events inside specific fiducial areas whose contours are 
defined by parameterizations of the kinematic variables 
of each particle. Within these well-defined regions, ac- 
ceptances and track reconstruction efficiencies are well 
understood using Monte Carlo simulations. 

After all selections have been applied, there remain 
about 130,000 exclusive pn + ir~ events. Fig. [S] shows the 
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FIG. 5: The distribution of squared missing mass 
M 2 {e'pTv + X) (GeV 2 ). The arrows show the exclusivity cut. 



V. CHARGED DOUBLE PION 
ELECTROPRODUCTION CROSS SECTIONS 



The kinematics of the 3-body pir + Tr~ final state are 
unambiguously determined by five independent vari- 
ables |64[. However, the choice of these variables is not 
unique. In this section we specify the kinematic variables 
we employ to describe the p7r + 7r _ final state and related 
phase space element for the 5-fold differential cross sec- 
tion. Note, that the double pion cross sections for vir- 
tual photon absorption are 5-fold differential, while the 
double pion electroproduction cross sections are 7-fold 
differential, since they contain the additional variables 
W and Q 2 . Then we describe the procedure to evaluate 
the 5-fold differential charged double pion cross section 
from experimental data from the 7-dimensional event dis- 
tributions. The 5-fold differential cross section contains 
complete information on double pion production at fixed 
Q 2 and W. However, the limited statistics do not al- 
low direct study of the 5-fold differential cross sections. 
Therefore, for the physics analysis we use various 1-fold 
differential cross sections, obtained by integrating the 5- 
fold differential cross section over the four other variables. 



Q 2 versus W distribution for the selected 2ir events. 
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FIG. 6: (color online) Q 2 (GeV 2 ) versus W (GeV) distribu- 
tion for the selected 2tt events. The grid shows the binning 
used for the evaluation of the cross section. Only cells inside 
the allowed phase space were used for that purpose. 
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FIG. 7: Kinematic variables for the ep — > e'p'ir + iv~ reaction 
(choice 2 in Sect. IV A"| . The top plot shows the n~ spherical 
angles 9 n - and <p v - , while the bottom plot shows the angle 
-pH-rr+p') °f the plane defined by the momenta of the final 
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p tv pair with respect to the plane comprised by the momenta 
of the initial proton p and tt~ . 
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A. Kinematic variables 

We adopt the following set of variables to describe the 
3-body final state: 

• invariant mass of the first pair of particles M12; 

• invariant mass of the second pair of particles M23; 

• the first particle solid angle fi; 

• the angle between two planes: one of them (plane 
A) is defined by the 3-momenta of the virtual pho- 
ton and the first hadron. and the second plane 
(plane B) is defined by the 3-momenta of the two 
other hadrons (see Fig. [7J. 

We use three different assignments for the first, second 
and third final state hadrons: 

• invariant mass of the pn + pair, invariant mass of 
the ir + w~ pair, the final proton spherical angles 
Op and ip p and the angle a^pp')^-*-) between the 
two planes: B, composed by the momenta of the 
7r + 7r~ pair, and A, composed by the momenta of 
the initial and final protons (choice 1) ; 

• invariant mass of the 7r + 7r _ pair, invariant mass of 
the p-K + pair, 7r _ spherical angles 0^- and (p^- and 
the angle a( p7r -)( p / 7r +) between the two planes: B, 
composed by the momenta of the final state proton 
p' and 7T + , and A, composed by the initial state 
proton p and tt~ (choice 2); 

• invariant mass of the pir + pair, invariant mass of 
the pK~ pair, ir + spherical angles 6 V + and tp n + and 
the angle a^+^pi^-) between the two planes: B, 
composed of the momenta of the final state proton 
p' and 7T~ and A, composed by the initial state 
proton p and tt + (choice 3). 

The 5-fold differential cross sections were obtained for 
all three sets of variables. The emission angles for the 
final particles in the second set of variables are shown in 
Fig. [T] For the other sets the emission angles arc defined 
in a similar way. In the physics analysis, described in 
Sect. IVI1 the second set of variables is used. These vari- 
ables are suitable for the description of charged double 
pion electroproduction through a 7r~A ++ intermediate 
state, which represents the main contributor of all iso- 
bar channels in the kinematic area covered by our data. 
The relations between the four momenta of the final state 
hadrons and the kinematic variables may be found in Ap- 
pendix A. 



of the first pair M\ 2 and the second pair M23 of the final 
state particles, solid angle for the first final state particle, 
and the angle on between the two planes A and B. The 
cross sections were determined only for those W cells, 
that arc fully inside the kincmatically allowed area. Spe- 
cial procedures were developed in order to evaluate the 
1-fold differential cross sections for the final state parti- 
cle invariant mass values near edges of the reaction phase 
space (see Sect. lVDIVEp . All frame-dependent variables 
and the cross sections were evaluated in the center-of- 
mass frame. 

For the second choice of kinematic variables, the 
7-fold differential cross sections dWdxy 1 ^^ > = 
dM pn +dM n + 7r -dCl n -dai pn -]ip/ n +] are given by the num- 
ber of pn + ir~ events AiV and efficiencies e in the 7- 
dimensional cells as: 

da 

dWdQ 2 d 5 r 2 = 

-A (1) 

e-e ch -R AWAQ 2 A 5 t 2 L y ' 

The number of events inside the 7-dimensional cells were 
corrected for contamination from the target walls, which 
was measured in separate runs with an empty target cell. 
The efficiency e in any 7-dimcnsional cell was determined 
in detailed Monte Carlo simulations. The inactive zones 
of CLAS and all cuts on phase space used in the event 
selection were included in the efficiency evaluation. The 
factor e c h accounts for the Cherenkov counter efficiency, 
which was determined separately. R accounts for radia- 
tive corrections. The integrated luminosity L was deter- 
mined from the total beam charge Q measured in the 
Faraday cup, combined with the information on target 
length and target density: 

q e M H 

where q e is the elementary charge, D t is the density of 
hydrogen (D t = 0.073 g/cm 3 ), l t is the length of the 
target (l t = 5 cm), Mh is the molar density of hydrogen 
(Mh = 1 g'/mol), and Na is Avogadro's number. The 
luminosity value was verified by reproducing elastic ep 
cross sections with the same data set. A comparison of 
the elastic ep scattering cross sections determined from 
our data with a parameterization of world data given in 
Ref. [H| showed agreement within better than 5 %. AW 
and AQ 2 are bins over W and Q 2 , and A 5 t 2 represents 
the element of hadronic 5-dimcnsional phase space for 
the second choice of kinematic variables: 



B. Evaluation of charged double pion cross sections 



The selected double ir events were collected in 7- 
dimensional cells, composed of W, Q 2 , invariant masses 



A 5 r 2 = AM pw +AM x+x -Acos(6» 7r -)Av? 7r -Aa (p7r - )(p , 7r + ) . 

(3) 

The 7-fold differential cross sections for the other two 
choices of kinematic variables listed in Sect. IV Al mav be 
obtained from eq. ([3]), by substituting the phase space 
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element A 5 T2 with: 

A 5 rx = AM JOT +AM 7r+7r -Acos((9p)A^pAa(pp')( 7r - 7r +); 



(4) 



A 5 r 3 = AM pw -AM p7r +Acos(6' 7r +)A^ 7r +Aa( p7r +)( f 



In the single photon exchange approximation, the 7- 
fold differential electron scattering cross section is related 
to the hadronic 5-fold differential cross section as 1661 1: 



da 



d 5 T 2 " + P 



d 4 T 7r + n = dM^+^-dfl^-da 



da 



dM^-p 



(pir )(p'tt"'") i 



(8) 



d 5 T 3 * P 



do 



dM p7r + dM v + v -dQ n -da( jm -)( p > v +) 



T v dWdQ 2 dM p7r + dM v +n- dCl n - da^-^i^ 
where T v is virtual photon flux, given by 



(5) 



a 



W(W 2 - M 2 ) 



II - e)Q 2 



(6) 



and a is the fine structure constant 
ergy, M } 

transverse polarization given by 



Eb is the beam en- 
p is the proton mass, and e is the virtual photon 



1 + 2 1 



Q 2 



tan 



(7) 



where v is the virtual photon energy and 6 e is the electron 
scattering angle in the laboratory frame. W, Q 2 and 8 e 
were taken at their respective bin centers. 

The limited statistics do not allow use of correlated 
multi-fold differential cross sections for physics analysis. 
In the physics analysis we instead used 1-fold differential 
cross sections, obtained after integration of the 5-fold dif- 
ferential cross sections over four kinematic variables. The 
number of 5-dimensional bins contributing to the indi- 
vidual bins of the 1-fold differential cross sections, range 
from 375 at W = 1.31 GeV to 1600 at W > 1.38 GeV. 
Summing up events in all 5-dimensional bins, reasonable 
statistical accuracy is achieved for the 1-fold differential 
cross sections (see Figs. [9|H3| ??.fT9 l l20 |) . We obtained 
in each (W, Q 2 ) bin covered by measurements a set of 
nine 1-fold differential cross sections, consisting of ir + n~ , 
pir + , and p7r~ mass distributions, 0i angular distribu- 
tions, as well as 3 distributions over angles q, (i=l,2,3), 
where the index i stands for the z-th set of kinematic 
variables, defined in the Sect. IV Al These 1-fold differen- 
tial cross sections represent the integrals from the 5-fold 
differential cross section over four variables as: 



da 



d 5 <T 

d^7 2 



d 4 T,„ 



7r+7r- 



tf^+Tr- = dM v - p dfl n -da^ pv -^ plw +y, 



d 4 T^~ p = dM 7V + 1l -dFlT T +da(pi r +)(p>T l -)', 



do 



d(—cos8{) 



d 5 o 
d^ 



d n; d n = dMn+n-dMn+pdipidoti 



In the actual cross section calculations the integrals in 
cq. ((SJ were substituted by the respective sums over the 
5-dimcnsional kinematic bins for the hadronic cross sec- 
tions. 

All cross sections represent independent 1-dimensional 
projections of the 5-fold differential cross sections. Any 
of the nine 1-fold differential cross sections provides inde- 
pendent information and cannot be computed from the 
8 remaining projections. 



C. Interpolation of 5- fold differential cross section 
into the CLAS detector areas of zero acceptance 

As discussed in Sect. IIVI the CLAS detector has areas 
with zero acceptance. The contributions from the 5-fold 
differential cross sections in such areas must be taken into 
account in order to obtain the integrated 1-fold differen- 
tial cross sections. We developed a special procedure to 
extend the 5-fold differential cross sections in the areas 
of CLAS with zero acceptance. 

We used the general <pi (£=1,2,3) dependence of the 5- 
fold differential cross sections fixing the other 4 kinematic 
variables: 



d 5 o 

— — = A + Bcos2o9i + Ccos<z>i + 
d b n 

B sin2<^i + C sin<y9i . 



(9) 



The first three terms are valid for any exclusive chan- 
nel and for any kind of particular reaction dynamics, 
being a consequence of rotational invariancc of the 
production amplitudes. The last two terms appear in 
the 5-fold differential cross sections for 3-body final 
states. After integration over the on angles, these two 
terms vanish as a consequence of parity conservation. 
The statistics in the populated bins is too small to allow 
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FIG. 8: Comparison of fully integrated charged double pion cross sections, obtained with (squares) and without (open circles) 
accounting for contributions from the CLAS areas of zero acceptance, as described in the text of Sect. IV Cl 
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FIG. 9: Mass distributions of p-K (top) and ir tv~ (middle), 
and 7r~ angular (bottom) distributions, obtained from differ- 
ent ways of interpolating the 5-fold differential cross sections 
into the CLAS areas of zero acceptance. Open squares corre- 
spond to 1-fold differential cross sections estimated without 
the contributions from the CLAS areas of zero acceptance. 
1-fold differential cross sections, obtained with the Rj coef- 
ficients in eq. (|10|l calculated within the framework of the 
JM03 and JM05 models, are shown by triangles and stars, 



evaluation of the A, B, C, B , and C coefficients from 
the data in the populated bins alone. Hence, we used 
both data and input from the models fit to the data 
to evaluate these coefficients. The coefficient ratios Rj 

= B/A, R 2 = C/ A, R 3 = B'/A, i? 4 = C/A ) 
were taken from phenomenological models for charged 
double pion electroproduction fit to our data. The 
coefficient A was determined from the data on the 5-fold 
differential cross sections in the populated 5-dimcnsional 
bins 



diM 



- d.Qdct, 



E 

Aip _ 



^(pTT~ )(p'7T+) 



-A<p 



A(2ir - Aip) - RiA J cos(2(p w -)d(p w - 

Aip 

—R2A J cos(<y9 7r - ) dtfir- - R3A J sin(2p„~) d<p n - (10) 

Aip Aip 

-RiA J sin(</? w -) d(p w - , 

Aip 

where the sum is running over the populated 5- 
dimensional bins, while the integrals are taken over the 
CLAS areas of zero acceptance Aip, The 5-fold dif- 
ferential cross sections in the CLAS areas of zero ac- 
ceptance were estimated from eq. © with coefficients 
A,B,C,B',C' calculated from eq. (fTU|) . 

In order to determine the ratios Rj within the frame- 
work of the phenomenological models, we propagated 
the 5-fold differential cross sections into the CLAS ar- 
eas of zero acceptance, using the JM03 model predic- 
tions for the shape of the 5-fold differential cross sections 
[H, [13, HH . The parameters of the JM03 model were de- 
termined from previous CLAS charged double pion data 
in the resonance region (2l| . In this way preliminary 
estimates for the 1-fold differential cross sections were 
obtained. Similar approaches to propagate the 5-fold 
differential cross sections into the CLAS of zero 

acceptance were used in previous charged double pion 
data analyses, published in Refs. [2l|, [24j]- In the next 
step, the parameters of JM03 were further adjusted to 
reproduce preliminary estimates of the 1-fold differential 
cross sections. The Rj coefficients were calculated within 
the framework of the JM03 approach after mentioned ad- 
justment of the JM03 parameters. The coefficients A for 
the if independent parts of the 5-fold differential cross 
sections (see eq. ((9])) were obtained from the data in the 
populated bins, according to eq. (JTUJ), using the improved 
estimates for Rj. Finally, the 1-fold differential charged 
double pion cross sections were obtained as described in 
Sect. IV Bl using the 5-fold differential cross sections in 
the CLAS areas of zero acceptance determined from eq. 
© with values of A, B, C, B and C determined as 
described above. 
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FIG. 10: (color on line) W (top panel) and Q 2 (bottom panel) 
distributions. The black lines represents data, and the red 
lines simulation. 



Since the model was used to evaluate the ratios Rj , the 
model assumptions affect mostly the (^-dependent parts 
in eq. ©■ After integration over tpi angles, these parts 
disappear. Nevertheless, as follows from eq. (flO|) . the 
model assumptions used to interpolate the charged dou- 
ble pion cross sections into the CLAS areas of zero ac- 
ceptance will increase the uncertainties of the 1-fold dif- 
ferential cross sections obtained in our analysis. Detailed 
studies to evaluate these uncertainties were carried out 
and are described below. 

First, we estimated the overall contribution from the 
inefficient areas to the 1-fold differential cross sections, 
calculating them in two ways: (1) by including the con- 
tributions from the areas of zero acceptance in CLAS, 
as described above, and (2) by excluding them. In all 
cases, the two sets of values were found to overlap well 
inside the statistical uncertainties in the entire kinematic 
region covered in the experiment. One example is shown 
in Fig [9] 

In Fig [5] we show a comparison of the fully integrated 
charged double pion cross sections. Again, within the en- 
tire kinematic area the differences between the two sets of 
cross sections are well within the statistical uncertainties. 

In the next step we investigate how the 1-fold differen- 
tial cross sections may be affected by the model assump- 
tions used in the procedure described above. Since only 
(^-independent parts of eq. © contribute to the 1-fold 



differential cross sections, we need to know the model 
uncertainties just for the A coefficients. As follows from 
eq. j9]), the influence of the ^-dependent parts on the 
A coefficients depend on: a) the relative contributions 
of the (^-dependent parts to the 5-fold differential cross 
sections; b) the ratio A</5/27r, where Aip is the overall 
(p coverage of the CLAS areas of zero acceptance. The 
relative contributions from the ^-dependent parts is es- 
timated by fitting the tfi angular distributions W , using 
eq. ((9]) with A, B and C coefficients as free parameters. 
The ipi angular distributions were obtained as integrals 
from the 5-fold differential cross sections over the other 
four variables. The B' and C terms should be equal to 
zero since these terms are integrated over the on angles. 
The contributions from B and C in eq. range from 10 
% to 50 %. For the majority of bins in Q 2 and W, these 
contributions range from 15 % to 25 %. The upper limit 
for the model dependence of the A coefficients has been 
estimated by replacing in eq. (|10|) all cosines by unity 
and by assuming A(p/2ir ^0.2 for the geometrical cover- 
age of the zero acceptance areas. With these assumptions 
we can calculate the model uncertainty of the A coeffi- 
cients, as the product of the maximal contribution from 
the ^-dependent parts to the 5-fold differential cross sec- 
tion (0.5) I 2 ' and the geometrical coverage of the CLAS 
areas of zero acceptance (0.2), resulting in an upper limit 
of 10 %. However, this limit was obtained with extremely 
conservative estimates for the integrands in eq. (|10p , the 
CLAS areas of zero acceptance Ap/2ir, and the relative 
contributions of the (^-dependent parts. More realistic 
estimates, outlined below, result in uncertainties of a few 
percent. 

In Fig. [9] we compare the results obtained using vari- 
ous models to estimate Rj . The JM03 and JM05 models 
are rather different in the description of the 5-fold dif- 
ferential cross sections. The JM05 approach provided a 
much improved treatment for the direct charged double 
pion electroproduction mechanisms |4lj |. It also con- 
tains an additional contact term that was introduced to 
improve the description of the ttA isobar channels. The 
interpolations of the 5-fold differential cross sections into 
the inefficient areas using these two models for Rj, and 
in addition taking off the contributions proportional to 
sin 2ip and sin tp, result in minor modifications well inside 
the statistical uncertainties. 

Finally we eliminate the contributions from the in- 
dependent parts in eq. (fTU)l and estimate the A coef- 
ficients from data in the populated bins. The results are 
shown in Fig. [9] (diamonds). Again, the estimated cross 
sections are well inside the statistical uncertainties of the 
data. 



[1] i=l,2,3 and stand for the set of kinematic variables, defined in 
Sect lVAl 

[2] ratio sum of (^-dependent parts over the full cross section. 
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FIG. 11: Comparison between measured (solid lines) and simulated (dashed lines) event distributions for various final state 
variables. (W = 1.4125 GeV, Q 2 = 0.525 GeV 2 ). 



D. Event reconstruction efficiencies 

A Monte Carlo event generator [63| was used to evalu- 
ate the event reconstruction efficiencies. The event gen- 
erator contains the main meson production channels in 
the resonance region. The efficiency for detection of the 
pir + iT~ in the final state was studied in detailed simula- 
tions that included the 2tt as well as 37r final states. The 
latter was needed to account for multi-pion background 
in the selection of charged double pion events, when ap- 
plying exclusivity cuts. These events were processed us- 
ing the same reconstruction program M , event selection 
procedures and fiducial cuts as for the events collected in 
the experiment. Efficiencies in the 7-fold differential bins 
were determined as the number of reconstructed events 
over the number of generated events, and used in eq. ([1]) 
to evaluate the 5-fold differential cross sections. 

As can be seen from Figs. [TU] and [TTJ the W and Q 2 
dependencies are well reproduced, as are major features 
in the event distributions over the final state hadronic 
variables. Differences between the measured and simu- 
lated distributions seen for the n + p invariant mass and 
7r~ angular distributions have little impact as efficiencies 
inside these areas are smooth. The event generator is 
therefore adequate to evaluate the event reconstruction 
efficiency for major parts of the kinematic range covered 



in the experiment. 

The limited phase space available for events in the low 
mass region with W < 1.40 GeV requires a different ap- 
proach for determining the event reconstruction efficien- 
cies for the invariant mass distributions. In this region we 
found rapid variations of efficiency from invariant masses 
inside the bins at the lowest edges for all mass distribu- 
tions. The use of an event generator that closely reflects 
the measured distributions is crucial in this area, where 
no 27r elcctroproduction data were previously available. 
Therefore, we have used an iterative procedure starting 
with the model event generator described above, and ex- 
tracted approximate cross sections for the different mass 
distributions. These were then used as a realistic in- 
put into the generated event distributions over invari- 
ant masses W gen (M k ) {k— r K + p 1 Tr~Tt + , n~p) for an ac- 
curate determination of event reconstruction efficiencies 
and cross sections. 

The improved estimates of event reconstruction effi- 
ciencies for the mass distributions ei mp (M k ) were ob- 
tained as: 



,(M fc ) 



W meas {M k )e{W,Q 2 
W gen (M k ) 



(11) 



[1] The correction factor that accounted for the events eliminated 
by the cut on the number of photoclcctrons in the Chcrcnkov 
counter was applied for the measured events only. 



where W meas (M k ) are event distributions in the invari- 
ant mass Mk and taken from the data. Both measured 
W m eas{M k ) and generated W gen (M k ) (k=ir+p event 
distributions were normalized to unity. The quantity 
e(W, Q 2 ) is the event reconstruction efficiency in a par- 
ticular (W, Q 2 )-bin, and was estimated using the event 
generator. 
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A comparison of the mass distribution obtained using 
the event generator with those estimated from eq. (fTTj) 
after further corrections, described in the Sect. IV E| is 
shown in Fig. [T3l 



E. Corrections for mass distributions 

After all of the previously discussed acceptance cor- 
rections, several mass distributions needed further cor- 
rections in order to account for the rapid variation of the 
cross sections inside some of the mass bins. For these 
bins, the cross sections were re-evaluated, using a bin- 
ning size reduced by a factor of 4. Cross sections at the 
nominal grid were compared to those obtained at the 
grid of reduced bin size and interpolated into the nomi- 
nal grid. In case of discrepancies, interpolated values of 
the cross sections were used, since they were determined 
with better mass resolution. 

A special procedure was developed to evaluate the 
cross sections at the smallest invariant masses using con- 
straints on the amplitude behavior near the phase space 
limits. All mass distributions at the smallest invariant 
masses were re-evaluated using a binning size reduced by 
a factor of 4 and interpolated over invariant masses in 
a way compatible with general requirements on a power 
low amplitude behavior near threshold: 



1 st bin 2 nd bin 



da 



o, 

C(M^- n+ -0.276) 

•o, 



dM k 

1.076)", 



lib 
"GeV 



M„- P (n+ P ) > 1.076 GeV 



M~- 



< 1.076 GeV 



> 0.276 GeV 



(12) 



M„+ w - < 0.276 GeV 



where C and a are free parameters fit to the cross sec- 
tions obtained with better mass resolution. The cor- 
rected cross sections were not evaluated in the mass areas 
closest to the threshold, which were affected considerably 
by the event migration. The bins of smallest invariant 
masses cover the mass intervals from the left edge of the 
next-to-smallest mass bin of reduced size to the right 
edge of the smallest mass bin of regular size (Fig. [T2")) . 
Differential cross sections in these bins were computed 
as: integrals from interpolating curves inside the bins di- 
vided by the bin size AM. The two solid curves in Fig.[T2l 
represent interpolating curves fit to the upper and lower 
boundaries of preliminary cross sections obtained with 
these improved mass binnings. Differences in the cor- 
rected cross sections, calculated using these two interpo- 
lations, give us the systematic uncertainties. These cor- 
rections only affect the lowest mass bins near the phase 
space limit. At larger invariant masses the cross sections 
were determined with the nominal bin size as described 
in the Sect. EE 

The comparison of mass distributions before and after 
all corrections described in Sect. IV Dl IV El is shown in 
Fig EH 
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FIG. 12: Corrections for mass distribution in the cross sec- 
tions near threshold. The filled circles represent the prelimi- 
nary cross sections obtained with the nominal bin size. The 
lowest mass bin of nominal size is shown by the solid vertical 
fines. Preliminary cross sections obtained using a bin size re- 
duced by a factor of 4 are shown by the open squares. The 
mass bins of reduced size are shown by the vertical dashed 
fines. The bin of lowest invariant masses AM, where cor- 
rected cross sections were evaluated, covers the interval from 
the dashed to the solid lines connected by the double sided ar- 
row. Corrected cross sections are given by the integrals from 
interpolating curves inside the AM bin over the bin width. 



F. Radiative corrections 

Radiative processes were evaluated using the proce- 
dure of Mo and Tsai [13] developed for inclusive pro- 
cesses, and incorporated in the event generator [63|. Ap- 
proaches that are capable at describing radiative pro- 
cesses in exclusive 2tt electroproduction, are not yet avail- 
able. 

The radiative correction factor R was determined as: 



R = 



N, 



norad 



(13) 



where N ra( i and N nora( i are the numbers of generated 
events in each (W, Q 2 ) bin with radiative effects switched 
on and off, respectively. This factor R was used in eq. (Q]) 
for calculations of the 5-fold differential cross sections. 

The particular hadronic tensor for exclusive 27r electro- 
production has impact mostly on radiation of the hard 
photons by the ingoing and scattered electrons. More- 
over its influence on observablcs decreases after integra- 
tion over the final state kinematic variables foa ]. The 
inclusive procedure for radiative corrections represents a 
reasonable approximation for the case of our data, since 
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W = 1.3375 GeV, Q 2 = 0.325 GeV 2 
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M,- p (GeV) 
W = 1.3625 GeV, Q 2 = 0.325 GeV 2 

T 




1.25 1.3 

M„- p (GeV) 

FIG. 13: Comparison of mass distribution cross sections be- 
fore (filled circles) and after (open squares) improvements, 
described in Sect. IVDIV El 



tions are reduced by factors of 2 to 4 after integration 
over the if angle for the emitted pion. Therefore, inte- 
gration over four variables in the case of charged double 
pion electroproduction is expected to reduce the radia- 
tive correction factor considerably, at least by a factor of 
4 |g| . 

Radiative corrections in the kinematics of this mea- 
surement were found to be less than 20 %. Therefore the 
contribution due to hard photon emission to the radiative 
corrections should be less than 10 %. They should be fur- 
ther reduced by more than a factor of 4 after integration 
of the 5-fold differential cross sections over four variables. 
Even a large uncertainty of 100 % in the contributions of 
hard photons would result in uncertainties of the overall 
radiative corrections to charged double pion cross section 
of only a few percent. 

The uncertainty in determining the cross sections 
caused by using the inclusive approximation for radia- 
tive corrections is well below the statistical uncertainties 
of the data, and is included in the systematic uncertain- 
ties 1 . 

W = 1.4625 GeV, Q 2 = 0.425 GeV 2 




1.2 1.4 

M n+p (GeV) 



0.4 0.6 

AW (GeV) 



we applied an exclusivity cut, which restricts the hard- 
ness for the emitted photons, and all 1-fold differential 
cross sections represent integrals of the 5-fold differential 
cross sections over four kinematic variables. 

We found that the relative contributions of hard pho- 
tons to R varied from 30 to 50 %. It is only this contri- 
bution that could be affected by the hadronic tensor and 
may be different in various exclusive channels. 

In Ref. [68| the effect of integration over kinematic vari- 
ables for the case of the exclusive single pion electropro- 
duction was studied, for which radiative processes have 
been evaluated exactly with the hadronic tensor derived 
from the fit to data. It was found that radiative correc- 



FIG. 14: Comparison between n + p and n~ir + mass distri- 
bution cross sections, obtained by integration of the 5-fold 
differential cross sections for various choices of the final state 
kinematic variables (choices from 1 to 3 are shown by file cir- 
cles, open squares and triangle respectively), defined in the 

Sect. Eg 



[1] The information obtained on the hadronic tensor from our JM 
model I4H based on a fit to the charged double pion data rep- 
resents valuable input for the future development of a fully ex- 
clusive radiative correction procedure for double pion electropro- 
duction. 



16 



TABLE I: Summary of systematic uncertainties for the fully 
integrated cross sections. The values represent averages over 
the kinematic areas covered by the data. 



Source 


Systematic 
uncertainty 


Integration 


3 % 


Fiducial cuts 


3 % 


Missing mass cut 


2.5 % 


Normalization 


< 5 % 


Binning effects 


3 % 


Radiative correction 


< 3.5 % 


Acceptances 


3 % 


Total 


< 9 % 



G. Results and systematic uncertainties 

In our analysis we determined nine 1-fold differential 
cross sections in each (W, Q 2 ) bin at invariant masses 
of the hadronic system from 1.30 to 1.57 GeV and at 
photon virtualitics from 0.2 to 0.6 GeV 2 with bins in W 
of 25 McV and in Q 2 of 0.05 GeV 2 . The data consist of 
7r + p, 7r~7r + , and ir~p invariant mass distributions, as well 
as 7r~, 7r + , and proton angular distributions, and 3 dis- 
tributions over angles aj (i=l,2,3) defined in Sect. IV A"l 
The full data set, consisting of 4695 cross sections, may 
be found in Ref. [23| ■ Fully integrated cross sections are 
shown in Fig. [T5] and Fig. 1161 These results represent 
the first comprehensive data set for charged double pion 
electroproduction at Q 2 < 0.6 GeV 2 . With respect to 
the previous data [2(j, the bin size in W is reduced by 
almost a factor of 10, while the binning in Q 2 is reduced 
by a factor of 5. 

Systematic uncertainties averaged over the kinematic 
range covered by the data are presented in Table [T] In 
the following we discuss the various contributions to the 
systematic uncertainties in more detail. 

To estimate uncertainties in the integration procedure, 
we compared the values of the fully integrated cross sec- 
tions obtained by integration over three sets of kinematic 
variables, defined in Sect. IV Al The ratios of r.m.s values 
for the integrated cross sections over their mean values 
were treated as systematics uncertainties related to the 
integration procedure. They are maximal at the lowest 
W value for all photon virtualities, and range from 5 to 
7 %. As W increases, they drop to ~ 1 % at the highest 
W value. 

In order to estimate the uncertainties in the integra- 
tion procedure for the 1-fold differential cross section, we 
compared their values determined in integration of the 
5-fold differential cross sections, obtained with three dif- 
ferent choices of the final state variables, described in 
Sect. IV Al When calculating these integrals for various 
kinematic variables, different sets of 5-dimensional bins 
contribute to the respective integrals. The efficiency was 
estimated for each set independently. Moreover, differ- 



ent inefficient areas contribute to the same cross sections 
estimated from integration over various kinematic grids. 
Therefore, a comparison of fully integrated and 1-fold dif- 
ferential cross sections, obtained by integration over var- 
ious kinematic variables, allows us to check the accuracy 
of the detector efficiency evaluations and propagation of 
the 5-fold differential cross sections in the CLAS areas of 
zero acceptance. 

Each of the three kinematic grids, discussed in 
Sect. IV Al contains the 7r + p invariant mass distribution, 
while just two grids contain 7r~7r + invariant masses. The 
angles describing the final state particles have unique 
assignments for each of three kinematic grids. We can 
therefore compare the results of integrations over 3 dif- 
ferent kinematic grids for the -n + p mass distributions. 
For the 7r _ 7r + mass distributions the integration over 
two grids can be compared. We found that these 1- 
fold differential cross sections coincide well within their 
statistical uncertainties and in the entire range of kine- 
matics covered by measurements. As an example, in 
Fig. ?? we show the comparison of mass distributions 
at W = 1.41 GeV and Q 2 = 0.425 GeV 2 , with the ones 
obtained from integrating the 5- fold differential cross sec- 
tions over different sets of kinematic variables. The com- 
parison of the fully integrated cross sections is shown in 
Fig. [15] The results differ by only a fraction of the sta- 
tistical uncertainties. 

The main contributions to the uncertainty in the over- 
all cross section normalization are given by uncertainties 
in the integrated luminosity and the electron detection 
and reconstruction efficiencies. These contributions have 
been estimated by measuring the well known elastic ep 
scattering cross sections. The comparison with a param- 
eterization [65| of the world data shows that the overall 
normalization is within a ~5 % uncertainty. 

In order to evaluate the systematics involved in defin- 
ing the final state exclusive process, we varied the missing 
mass cut used to identify the unmeasured tt~ , and mod- 
ified the fiducial regions where final state particles are 
selected. The average uncertainties are shown in Table fl] 

In Sect. IV CI we concluded that the contributions of 
the zero acceptance regions in CLAS affect the extracted 
differential and integrated cross sections well within the 
statistical uncertainties. Systematics uncertainties re- 
lated to propagating the data into the inactive areas of 
CLAS were estimated, assuming a 50 % uncertainty in 
the extrapolation of the 5-fold differential cross sections, 
resulting in 2 % to 5 % uncertainties, and increasing in 
Q 2 . 

The global systematics for radiative corrections, listed 
in the Table HI were calculated assuming the individual 
contributions are uncorrelatcd. The factor R, obtained 
from our Monte Carlo simulation (Sect. IV F[) . revealed 
no Q 2 dependence in the entire kinematic range of our 
measurements. The root-mean-square values for the R 
factors calculated at various Q 2 were assigned to the un- 
certainties for the radiative correction factor. Based on 
the estimates described in Sect. IV Fl we assigned an up- 
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FIG. 15: Comparison between fully integrated charged double pion cross sections, obtained in the integration of the 5-fold 
differential cross sections over three various choices of the final state variables, described in the Sect. IV Al The choices for 
kinematical variables from 1 to three are shown bv triangles, the open and full squares respectively 
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FIG. 16: Fully integrated 2n cross section at various Q 2 . Crossed areas represent systematic uncertainties. Full calculations 
within the framework of the JM06 model [4^] are shown by solid curves. The contributions from s-channel resonances and from 
non-resonant mechanisms are shown bv the dot-dashed and dashed curves respectively. 
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per limit of 2.5 % to the uncertainties related to the hard 
photon emission. 

The overall systematic uncertainty shown in Table U 
was obtained as the square root of the quadratic sum 
over the individual contributions. Applying the described 
procedures in each bin of W and Q 2 individually, we ob- 
tain the systematic uncertainties for the integrated cross 
sections as shown by the shaded areas in Fig. [16] 



VI. PHENOMENOLOGICAL ANALYSIS 

The comprehensive information on the 1-fold differ- 
ential cross sections for charged double pion production 
provided by the data enabled us to carry out a com- 
bined analysis of the nine 1-fold differential cross sec- 
tions in each W and Q 2 bin, and to establish all signif- 
icant mechanisms contributing to charged double pion 
electroproduction in the range W < 1.6 GeV and Q 2 = 
0.2-0.6 GeV 2 . Before our data were available, the infor- 
mation on 2tt electroproduction mechanisms in this kine- 
matic area was rather scarce and quite uncertain. Rea- 
sonable knowledge on 2ir electroproduction mechanisms, 
achieved in our data analysis, also allowed the separa- 
tion of resonant and non-resonant contributions needed 
for the evaluation of nucleon resonance electrocouplings. 

The presence and strengths of the contributing 2tt elec- 
troproduction mechanisms were established by studying 
the kincmatical dependencies in differential cross sections 
and their correlations in a variety of available obscrv- 
ables. The analysis was carried out using a phenomcno- 
logical model developed in the past few years by the Jef- 
ferson Laboratory - Moscow State University ( JM) col- 
laboration [H, Ell, [H, 0, [H. Within the JM model 
developed up to 2005, called JM05, the major part of 
7r + 7r~ production at W < 1.6 GeV is due to the contri- 
bution from the ttA isobar channels. This is supported 
by the data in that the peaks from A ++ (1232) reso- 
nance were clearly seen in all n + p mass distributions at 
W > 1.4 GeV, while other mass distributions did not 
show any visible structures. The contributions from all 
other isobar channels included in JM05 become negligi- 
ble in this deeply sub-threshold region of W < 1.6 GeV. 
Analysis of the recent CLAS 2ir data, presented in this 
paper, allowed us to study manifestations of the JM05 
mechanisms and search for other contributing processes 
in the still unexplored kinematics area of photon virtual- 
ities from 0.25 to 0.6 GeV 2 . This analysis also enabled us 
to access the dynamics of the contributing mechanisms 
at a phenomcnological level, fitting them to all measured 
observables combined. Our knowledge on the contribut- 
ing mechanisms was extended considerably, resulting in 
the recent version of the JLAB-MSU model, which we 
refer for as JM06. A detailed description of the JM06 
model version may be found in a separate paper [44| . 
Here we discuss the basic ingredients of JM06 and the 
major results. 

The 7*p — > 7r + 7r~p production amplitude within the 



JM05/JM06 model versions can be written as: 

rp th7tA i rpdir /-i a\ 

J-yN,inrN — ± 1 * N ^nN + 1 1 *N,-k^N 

with 

^7*Ar,7r7rAf = [*7»JV,ttA + *t*A>A + ^7*JV,tta] ' (15) 

where Ga is a propagator of the A intermediate state, 
and the vertex function Ta.kN describes the A(1232) — > 
ttN decay. The mechanisms in the above equations for 
the JM05 model version are illustrated in Fig[T7l The di- 
agram — > N* , A* — » 7rA in the second row of the fig- 
ure is the resonant term t^, N ^ m ec L- (©■ It is param- 
eterized as a Breit-Wigner form (3|| and calculated from 
all well established N* states with masses less than 2.0 
GeV that have hadronic decays to the irirN final states. 
In the kinematic area covered in our measurements, only 
the Pn(1440) and Di3(1520) nucleon resonances have 
the contributions in 1-fold differential cross-sections, that 
are outside of the data uncertainties. The non-resonant 
term t^°JJ L 7rA is calculated from the well established Born 
terms of 7* N — > 7rA [2^, [36[. Their amplitudes are pre- 
sented in Ref. [36|, [44[ . The final state interactions were 
treated effectively within the framework of the absorp- 
tive approximation [36| . Additional contact terms were 
implemented in order to account for possible contribu- 
tions from other mechanisms to 7rA production, as well 
as for hadronic interactions of 7rA states with other open 
channels [4l|, [7(| • These extra contributions to the ir A 
isobar channels were previously established [4l| in anal- 
ysis of CLAS 2tt data [2l| and confirmed by the data 
of this paper. A parametrization for these amplitudes 
will be presented in Ref. [HI ■ The contributions from all 
isobar channels combined account for from 70 to 90 % 
of the 2ir fully integrated cross sections in the kinematic 
area covered in our measurements. A remaining part 
of cross sections comes from direct charged double pion 
production mechanisms, when the p7r + 7r~ final state is 
created without the formation of the intermediate quasi- 
two-body states with unstable hadrons. In the JM03 
model version we started with a simplest parametriza- 
tion for direct 2-7T mechanisms as the three body phase 
space with the amplitudes that were independent from 
the final state kinematic variables and fit to the data in 
each bin of W and Q 2 independently [3?], H8|. However, 
this parametrization did not allow us to reproduce steep 
dependencies in ir~ angular distributions at the back- 
ward angles, clearly seen both in the previous CLAS 2ir 
electroproduction data [2l[ and in the data of this paper. 
The example is shown in Fig. 1191201 In order to reproduce 
such data behavior, the direct 27r production mechanisms 
were parameterized in JM05 [4fj| . [43[ in terms of a contact 
vertex and a unspecified particle-exchange amplitude, de- 
scribed by the effective propagator, which depends expo- 
nentially from a running four-momenta squared. The di- 
agrams in the bottom of Fig[T7]represents the direct term 



20 






jc + (7C-)(jc-) y 

p'(TC + )(p') 

+ 

Jt-(p')(3I + ) P 

JM05 parametrization 



p'(3I + )(7f) 




FIG. 17: The mechanisms of JM05 model [4Cj, |4l| contributing to 2-ir electroproduction at low W and Q 
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FIG. 18: Direct charged double pion production mechanisms 
in JM06 model. 



T^l r N ^ N in eq-CEU) that was introduced in Refs.galll 
to describe the direct 7^V — > irirN mechanisms. The pa- 
rameterization of this term in JM05 was given explicitly 
in Ref. gO]. 

This parametrization allowed us to describe success- 
fully the previous CLAS 2-7T data [2l[ and the data of 



our paper on three invariant masses and 7r~ CM. an- 
gular distributions. Analysis of nine 1-fold differential 
cross-sections, which was carried out for the first time, 
demonstrated the shortcomings in description of n + and 
p CM. angular distributions related to the parametriza- 
tion of the direct 2ir electroproduction in the JM05 model 
pT0j | . Our data offer a compelling evidence for neces- 
sity to modify a description of the direct 2tt produc- 
tion mechanisms. In the recent JM06 model version 
they were parameterized as two subsequent unspecified 
particle-exchanges amplitudes, shown in Fig. [TS] for var- 
ious assignments of the final state particles. Both prop- 
agators describing unspecified particle exchanges were 
parametrized by the same exponential dependence from 
a running four-momenta squared. Explicit parametriza- 
tion for the amplitudes of these mechanisms will be pre- 
sented in a separate paper of Ref. [44[ . 

Within the framework of the JM06 approach we were 
able to describe the 27r data of our paper in the en- 
tire kinematic range covered by the measurements. As 
a typical example, the model description of the nine 1- 
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FIG. 19: Description of CLAS charged double pion differential cross sections at W = 1.51 GeV and Q 2 = 0.425 GeV 2 within the 
framework of the JM06 model. The full calculations are shown by the solid lines. The contributions from 7r~A ++ and 7r + A° 
isobar channels are shown by the dashed and dotted lines, respectively. The contributions from direct charged double pion 
production processes are shown by the dot-dashed lines, a; angular distributions were calculated with the JM06 parameters 
fit to the other 6 differential cross sections. 



fold differential charged double pion cross sections at 
W = 1.51 GeV and Q 2 = 0.425 GeV 2 is shown in Fig.[]J 
together with the contributions of all mechanisms incor- 
porated in the JM06 description. 

The shapes of the cross sections for the different mech- 
anisms are substantially different in the various observ- 
ables, but highly correlated by the reaction dynamics. 
Moreover, we found no need to implement additional 
mechanisms beyond the ones already included in JM06. 
Therefore, the successful description of all nine 1-fold dif- 
ferential charged double pion cross sections allowed us to 
identify all significant contributing processes and access 
their dynamics at the phenomcnological level. To check 



the robustness of the results obtained within the frame- 
work of the JM06 model, we fit the model parameters to 
a limited set of data that included only six differential 
cross sections: all three invariant masses and three an- 
gular 9i (i=tt~ , 7r + , p) distributions for the final state 
hadrons. The remaining three distributions over the 
ai angles were calculated, keeping the JM06 parameters 
fixed. A good description of all of the ai distributions 
was achieved throughout the kinematics covered by the 
measurements, giving us confidence that the main pro- 
cesses are described within the JM06 model. 

The separated resonant and non-resonant contribu- 
tions to the fully integrated cross sections are also pre- 
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FIG. 20: Resonant (dot-dashed lines) and non-resonant (dashed lines) contributions to the charged double pion differential 
cross sections at W=1.51 GeV and Q 2 =0.425 GeV 2 . The full JM06 calculations are shown by black solid lines. 



sented in Fig. 1161 The resonant and non-resonant parts 
of the differential cross sections at W = 1.51 GeV and Q 2 
= 0.425 GeV 2 are shown in Fig. [20] Non-resonant mech- 
anisms represent a major contributor in entire kinematic 
area covered by our measurements. The relative resonant 
contributions increase with Q 2 but remains below than 30 
%. However, the resonant contributions increase rapidly 
at W > 1.4 GeV, where they become larger than the 
data uncertainties and may be determined from the data 
fit . The dominant part of the resonant amplitudes comes 
from the P u (1440) and D 13 (1520) states combined. The 
phase space limitations prevent P 33 (1232) decays to the 
final states with two pions. There is also no evidence 



for substantial decays of the 5 l n(1535) resonance with 
two pion emission, while the tail from the nucleon exci- 
tations with masses above 1.6 GeV is well inside the data 
uncertainties. 

The shapes of the resonant and non-resonant contri- 
butions are quite different for most differential cross sec- 
tions, especially the angular distributions. Moreover, 
the correlations between kinematical dependencies of the 
resonant /non- resonant parts in various 1-fold differential 
cross-sections arc also quite different. Therefore, a rea- 
sonable description of all differential cross sections com- 
bined, which is achieved within the framework of the 
JM06 model, provides an evidence that a proper isola- 
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tion of the resonance contributions can be obtained from 
global fits to all available cross section data, even in a 
case of relatively small (< 30 %) resonant contributions 
to the fully integrated cross-sections. It may be seen in 
Fig. [201 where derived from the data fit the resonant con- 
tributions are presented. 

While the evaluation of N* electrocouplings is outside 
the scope for this paper, it is nevertheless noteworthy 
that our data open up the possibility to determine elec- 
trocouplings for Pu(1440) and Z?i 3 (1520) states at Q 2 
< 0.6 GeV 2 from double pion electroproduction for the 
first time. The preliminary results may be found in the 
Ref. [70l ]. The kinematical region of low photon virtu- 
alities is expected to be particularly sensitive to meson- 
baryon dressing effects in the resonance structure (45ll49j|. 

Moreover, in this kinematic region there are also CLAS 
data on single pion electroproduction [1, [1,0, SB III, III]- 
At W < 1.5 GeV only single and double pion exclusive 
channels contribute to the total meson production cross 
sections off protons. The amplitudes for various non- 
resonant mechanisms contributing to single and double 
pion electroproduction, obtained in a phenomenological 
analysis of CLAS data with the framework of the JM06 
model [HI, represent valuable information for the study 
of nucleon resonance transitions in combined analyses 
of these major channels. Advanced coupled-channel ap- 
proaches arc under development at the Excited Baryon 
Analysis Center (EBAC) at Jefferson Lab (MUl]. Phe- 
nomenological analysis of the double pion data of this 
paper, described in detail in Ref. [44j, is of particular 
interest for this activity. 



VII. CONCLUSIONS 

In this paper we have presented a large body of electro- 
production data for the process j v p — > pir + TT~ . The large 
acceptance of CLAS allowed extraction of the 1-fold dif- 
ferential and fully integrated charged double pion cross 
sections. The results were tested for robustness by using 
different integration grids, that showed that consistent 
results are obtained independent of the specific integra- 
tion procedure. 

1-fold differential and fully integrated cross sections 
were obtained for W from 1.3 GeV to 1.6 GeV and Q 2 
from 0.2 GeV 2 to 0.6 GeV 2 . The high statistics and good 
momentum resolution of the measurements allowed us to 
use bin sizes of AW = 25 MeV and AQ 2 = 0.05 GeV 2 , 
which are at least a factor of 5 smaller than the ones 
used in previous measurements. For the first time, nine 
independent differential cross sections in each bin of W 
and Q 2 were measured. 

The phenomenological analysis of the cross sections 
within the framework of the JM06 approach [43l . l44j al- 
lowed us to establish significant mechanisms contributing 
to charged double pion electroproduction for the kine- 
matics covered by our measurement. All differential and 
fully integrated cross sections obtained in our measure- 



ment can be reasonably described by the contributions 
from 7rA isobar channels and direct double pion pro- 
duction mechanisms, established from phenomenological 
analysis of our data. The analysis of resonance contri- 
butions to the differential cross sections indicates that 
these data are sensitive to the resonant parts of cross- 
sections, allowing us to study the excitation of Pi i (1440) 
and D 13 (1520) states by virtual photons off protons at 
Q 2 < 0.6 GeV 2 . 

The data also allowed us to determine cross sections 
and amplitudes for isobar channels, offering valuable in- 
formation for nucleon resonance studies in a global anal- 
ysis of the major meson electroproduction exclusive re- 
actions within the framework of an advanced coupled- 
channel approach currently under development in EBAC 
at Jefferson Lab [H, HI). 
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APPENDIX A: KINEMATICS VARIABLES FOR 
5-DIFFERENTIAL 2-tt PRODUCTION 
CROSS-SECTIONS 

In this appendix we present the final state kinematics 
for the second choice of variables defined in Sect. IV Al 
Since all momenta are measured in the lab frame, first 
we boost the 3-momenta of the final state particles in the 
cm. frame. All 3-momcnta used below, if not specified 
otherwise, are defined in the cm. frame. 



, Af T +„ and M w - V invariant masses were related 



to the four momenta of the final particles as: 



M w+7r - = y/(P„+ +P„-] 



AL 



7T + p' 



P. 



p p >Y 



(16) 



ir + , p) stand for the final state particle 



where Pi (i=n~ . 
four-momenta. 

The angle 9^- between the 3-momentum of the initial 
photon and the final state 7r _ in the cm. frame was 
calculated as: 



IP,- I IP* 



(17) 



The angle ip^ 
arctg 



was determined as: 



p 



arctg 



arctg 



arctg 



P, 



P _ 
P _ 

p _ 

1 T.7T 



P xn - > 0; P y „- > (18) 
f P xn - > 0; P^- < (19) 

(20) 



P, 



< 0; P m - < 



tt/2; P m - 

3tt/2; P m 



0; P 



J/7T- 

o- P 



P xw - < 0; 

> 
- < 0. 



P y7T - > 



(21) 



The calculation of the angle QY„— p )(„ 



(22) 
(23) 

between the 



two planes A and B (see Fig. [7j is more complicated. 
First we determine two auxiliary vectors 7 and j3. The 
vector 7 is the unit vector perpendicular to the P v - 
3-momentum, directed toward the vector —ri z and sit- 
uated in the plane composed by the virtual photon 3- 
momcntum and the 7r _ 3-momcntum P„- (see Fig. [7|)- 
il z is the unit vector directed along the z-axis (see Fig.[7|). 
The vector (3 is the unit vector perpendicular to the 3- 
momentum of tt~ , directed toward the tt + 3-momentum 
Pk+ and situated in the plane composed by the ir + and 
p' 3-momenta. Note that the 3-momcnta of the tt + , tt~ 
and p' arc in the same plane, since in the ccntcr-of-mass 
their total 3-momentum should be equal to zero. Then 
the angle between the two planes oi^-p)^^) i s: 



with the acos function running between zero and 7r, and 
the angle between the planes A and B running from zero 
to 2ir. To determine a in a range between n and 2-7T, 
we look at the relative orientation of the vector P v - and 
vector product S for the auxiliary vectors 7 and (3: 



7 x (3. 



If S is collincar to P_- 



(25) 



is determined from 



eq. (|24|) . In the case of anti-collinear vectors 8 and P^-: 
a^- p )(^+p') =2ir- acos(7^). (26) 

The vector 7 may be expressed through the particle 3- 
momenta as: 



7 = a a {-n z ) + b a rip _ 



1 



Vl-(nVK)) 2 
K = ~{np _ {-n z ))a a , 



(27) 



where rip _ is the unit vector directed along the 7r~ 3- 
momentum (see Fig. [7]). Taking scalar products ("/ftp _ ) 
and (77), it is straightforward to verify that 7 is the unit 
vector perpendicular to P n - . 

The vector f3 may be obtained as: 



[3 = aprip^ + bprip^ 



bp = -(n P ^ + rip^_ )ap , 



(28) 



where rip , is the unit vector directed along the 7r + 3- 
momentum. Again, taking scalar products ([3rip _) and 

(f3/3), it is straightforward to see that (3 is the unit vec- 
tor perpendicular to the 7r~ 3-momcntum. The angle 
a (ir-p)(ir+p') coincides with the angle between the vectors 
7 and (3. So, the scalar product (7/?) allows determina- 
tion of the angle a^-p^^+pt) in eq. (j2"4")) . 

The kinematic variables for other hadron assignments 
for the first, second and third final state particle de- 
scribed above, were evaluated in a similar way. 



C 71- p)( t < 



hp/) = acos (7/3), 



(24) 



